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NOMENCLATURE 

;mbol  Def init  ion  Dimensions 

A       Empirical  constant 

a       Empirical  constant  volts/  F 
B       Empirical  constant  ^       ^ 

,  Volt  S-.  ,      f  t  x  2 

C  Defined    by   Equation    17  ohms  sec 

C  Specific   heat    of   air  BTU/lb    °F 

P 

D  Diameter    of    hot-wire  inch 

E  Instantaneous   hot-wire/thermocouple    voltage  volts 

E  D.C.    component    of    hot-wire/thermocouple    voltage         volts 

E  Thermocouple    reference   voltage  volts 

e'  Fluctuating   component    of   hot 

wire/ thermocouple    voltage  millivolts 

h  Convective   heat    transfer    coefficient  BTU/hr'ft     •    F 

k  Thermal    conductivity    of    air  BTU/hr  •  f t  •  CF 

k  Time    average    thermal    conductivity    of    air  BTU/hr'ft'    F 


or 


k'       Fluctuating  component  of  thermal 

conductivity  of  air  BTU/hr ■ f t • ^F 

Z  Length  of  hot  wire  inch 

m       Empirical  constant 

Nu       Nusselt  number 

AP      Pressure  difference  between  points  inches  of 

upstream  and  downstream  of  orifice  plate  manometer  fluid 

Pr       Prandtl  Number  of  air 

2 
q"      Heat  flux  from  surface  of  hot-wire  BTU/hr -ft 


IX 


R       Pipe  radius  inch 

r        Distance  from  pipe  cenlerline  inch 

R       Hot-wire  resistance  at  local  air 

temperature  ohm 

R       Hot-wire  resistance  at  time  average 

local  air  temperature  ohm 

R       Hot-wire  reference  resistance  ohm 

o 

R       Hot-wire  resistance  at  operating 

temperature  ohm 


r '      Fluctuating  hot-wire  resistance  corresponding 
to  fluctuating  component  of  local  air 
temperature  ohm 

Re      Reynolds  Number 

S       Hot-wire  sensitivity  to  velocity 


u 


fluctuations  mi llivolts/f t /sec 


S  Hot-wire    sensitivity    to   temperature 
f luctuat  ions 

T  ,T  Instantaneous  local  air  temperature 
a 

T  ,f  Time  average  local  air  temperature 

3. 

T  Time  average  temperature  at  pipe  center  Line 

T  Wall  temperature 
o 

T  Wire  temperature 
w 

l  '  Fluctuating  component  of  local  air  temperature 

u  Instantaneous  local  flow  velocity 

u  Time  average  local  flow  velocity 

ii  Fluctuating  component  of  local  flow  velocity 

a  Hot-wire  resistivity  coefficient 

p  Instantaneous  air  density 

P  Time  average  air  density 

p  Fluctuating  component  of  air  density 


mi llivo 

It: 

5/°F 

°F 

°F 

°F 

°F 

°F 

°F 

ft /sec 

ft /sec 

ft /sec 

oF-l 

lb/ft3 

lb/ft3 

lb/ft3 

ABSTRACT 

Wiggins,  Edwin  George,  M.S.  in  Nuclear  Engineering,  Purdue 
University , January  1969.  An  Investigation  of  the  Use  of  the  Hot-Wire 
Anemometer  in  Non- Isothermal  Air  Flow.   Major  Professor:   Alexander 
Sesonske . 

Procedures  were  developed  for  the  calibration  of  the  hot-wire 
anemometer  in  a  non- isothermal  flow.   Various  methods  of  using  the  re- 
sults of  calibration  to  compute  turbulent  fluctuations  of  velocity  and 
temperature  were  considered,  and  the  Arya  and  Plate  modification  of  the 
Kovasznay  Fluctuation  Diagram  method  was  found  to  give  the  most  accurate 
results.   The  method  requires  independent  measurement  of  the  temperature 
fluctuations  which  was  obtained  with  a  rapid  response  thermocouple. 

The  calculated  values  of  velocity  turbulence  intensity  show  devi- 
ations of  36  percent  from  the  results  of  Laufer ;  however,  this  deviation 
is  to  be  expected  in  view  of  the  estimated  experimental  error.   These 
errors  were  probably  caused  by  contamination  of  the  wire  during  the 
lengthy  calibration  process. 

All  measurements  were  made  in  air  flowing  in  a  1.0625  inch  diameter 
pipe  at  Reynolds  Numbers  between  30,000  and  50,000.   A  heat  input  of  430 
watts  over  the  last  58  inches  before  the  probe  was  used  to  generate 
temperature  differences  between  wall  and  centerline  of  48  to  63°F  at  the 
probe  location. 

The  procedures  developed  and  the  experimental  problems  identified 
are  the  contribution  of  this  investigation  to  a  comprehensive  research 
program  of  non- isothermal  hot-wire  measurements  in  fluids  of  a  wide  range 
of  Prandtl  Numbers. 


INTRODUCTION 
The  hot-wire  anemometer  has  been  used  extensively  for  the  measure- 
ment of  velocity  fluctuations  in  isothermal  flow  of  gases.   The  principles 

12  3 

and  procedures  involved  are  discussed  by  Kovasznay    and  by  Kronauer  . 

4 
Corrsin  has  shown  that  in  addition,  the  hot-wire  anemometer  has  the 

potential  to  measure  simultaneous  velocity  and  temperature  fluctuations 

in  non- isothermal  flow.   Experimental  work  of  this  nature  has  been  done 

by  Corrsin  and  Uberoi  ,  Gibson  Chen  and  Lin  ,  Deissler  ,  Kunstman  ,  and 

9 
Arya  and  Plate  .   The  present  work  was  concerned  with  the  development  of 

experimental  procedures  and  the  identification  of  problems  associated  with 

making  of  such  measurements  in  a  non- isothermal  air  flow. 

In  isothermal  flow  the  only  flow  variable  affecting  the  anemometer 
output  is  the  velocity.   Thus  the  anemometer  voltage  fluctuation  is 
directly  proportional  to  the  velocity  fluctuation  and  the  anemometer 
steady  voltage  is  directly  proportional  to  the  mean  velocity. 

In  non- isothermal  flow,  however,  both  the  velocity  and  ambient 
temperature  fluctuations  affjpct  the  anemometer  output.   Considerable  dif- 
ficulty has  been  encountered  by  previous  workers  in  separating  the  two 

Q 

effects.   Errors  up  to  750  percent   have  been  reported  in  non- isothermal 

measurements,  while  in  isothermal  measurements  10  percent  error  or  less 
is  expected. 

The  non- isothermal  flow  measurements  of  the  present  work  were  pre- 

ceeded  by  velocity  fluctuation  measurements  in  isothermal  flow  with  a 


hot-wire  paralleling  the  work  of  Laufer    and  Sandborn  "  and  temperature 
fluctuation  measurements  in  non- isothermal  flow  with  a  rapid  response 
thermocouple.   Thus  the  velocity  and  temperature  fluctuation  intensities 
measured  with  the  hot-wire  in  non- isothermal  flow  could  be  compared  with 
results  of  more  direct  measurements. 


THEORETICAL  BACKGROUND 

The  constant  temperature  hot-wire  anemometer  consists  of  a  fine 
tungsten  or  platinum  wire  maintained  by  a  feedback  system  at  some  constant 
temperature  above  that  of  the  fluid  in  which  it  is  inserted.   The  con- 
vective  heat  loss  from  the  wire  is  balanced  by  joulean  heating  of  the 

wire  by  the  current  supplied  by  the  feedback  system. 

12 
King    solved  the  convection  problem  by  assuming  potential  flow 

around  the  wire  and  obtained  the  following  relationship: 

T  -T   m 
Nu  =  (A  +  BRe2)  (1  +  \     ™      &      )  , 

a 

where  Re  =   Pup  . 

Although  it  is  recognized  that  King's  assumptions  are  unacceptable,  the 
functional  form  of  his  relationship  has  proved  to  be  a  satisfactory  re- 
presentation of  experimental  data.   Thus  it  is  seen  that  the  joulean 
heating  current  (or  voltage)  required  to  maintain  the  wire  at  a  constant 
temperature  is  a  function  of  the  instantaneous  velocity  of  the  flow  and 
the  difference  in  temperature  between  wire  and  fluid. 

In  isothermal  flow,  the  temperature  difference  is  constant  and  the 
system  can  be  used  to  measure  the  mean  and  fluctuating  components  of  the 
instantaneous  velocity.   In  non- isothermal  flow  the  fluid  temperature 
has  a  fluctuating  component  that  affects  the  heat  transfer  in  addition 
to  the  effect  of  the  fluctuating  component  of  the  velocity.   Thus  in 


principle,    the    system  can   be   used    to  measure   velocity    and    temperature 
fluctuations,    since    it    is    sensitive    to  both    fluctuations. 


Hot-Wire    in    Isothermal   Flow 
The    following      parallels   Corrsin's        treatment    for    a   constant    current 
hot-wire    anemometer.      For   heat    transfer    to   flow    over    a   cylinder: 

1/9  1    T      -   T        m 

Nu   =    (A  +   B    Re^Ml   +  J  T  ]     •  O) 


In   the  range    of    interest  , 

,  T      -   T 

1  _w a      .  (2) 

2  T 


Thus  Equation  1  may  be  simplified  to 


1/2 
Nu  =  A  +  B  Re    .  (3) 


iy   definition, 


Nu  =  ™  = a^2 (4) 

k    k(T  -  T  ) 

w    a 

For  electrical  resistance  heating, 

,"  -  ^  • 

M  R   ttD£ 

w 

Substitution   of  Equations   4    and    5    into   3    yields 
2 


1/2 
(7Tk£A)    +    (iktB)    Re    ' 


By   definition, 


(5) 


R    (T      -   T    )  (6) 

ww  a 


Re   =  ^-   .  (7) 


Thus  , 


f2  n      1/2         1/9 

=    Uk£A)    +    (irk£B    [    ^   ]         )    u   ' 


R  (T   -  T  )  u 

w   w     a  \o ) 

The  quantities  in  parentheses  on  the  right  side  of  Equation  8  are 

constants.   For  constant  temperature  operation  of  the  hot-wire  in  an 

isothermal  flow  the  following  mean  and  deviation  quantities  are  used: 


u  =  u  +  u '     and     E  =  E  +  e ' . 

Substitution  of  Equations  9  into  Equation  8  yields,  after  small 
terms  have  been  discarded, 


72    i7    i  n   !/2   i  In  •   !/2 

R  /T  _T  n  =  (ttUA)  +  (nkV.B  [  *-  ]    )u    (i  +  —  ) 

w  w  a 


By  definition, 


T   -  T   =  -—    (R   -  R  ). 
w     a    aR    w     a 


According  to  the  binomial  theorem 


(9) 


(10) 


(11) 


,1/2       ,   , 

(1+*-)     =  l+±^  +  ....  (]2) 

u  Z    u  (12 ) 

Since  —      1,  it  is  assumed  that  higher  order  terms  may  be 

neglected  in  Equation  12.   Thus  Equation  10  becomes 

E  2Ee'  ,TTk£    ..,,Trk*    _  r£D,  1/",   -1/2 ,         u'      x 

R    (R-R   )+R    (R-R   )    =(^R-A)  +  (^"B[u]         >  (u        +  7T7I  }  '  (13) 

wwawwao  o  2u 


Time    averaging    of   Equation    13    yields 

E2  ,    Tik£  TikP    _    .    pD    ,1/2-l/2  „.. 

R   (R  -R   )    =    (  7^"  A)    +    (  ^R"  B    [   T   ]         )u        •  (14) 

w     w      a  o  o 


Subtraction   of  Equation    14    from    13    gives 


2Ee' 


R    (R      -    R    ) 

ww  a 


1    r  n        1/2 

aRo  M  2u1/2 


(15) 


or 


4Ee' 


—1/2 
R    (R      -   R   )    Cu   ' 

ww  a 


(16) 


Note    that 


c   =  JlMB    [   ££ 
uR  u 


1/2 


and    that    this    same   quantity    appears    as    the   coefficient    of   u      in 
Equation    14.      Squaring,    time    averaging,    and    taking    the    square 
root    of   Equation    16    yields 


4E/    e 


,2 


/   u'2 


-1/2 
R    (R      -    R    )    Cu    ' 

ww  a 


(17) 


(18) 


Rapid   Response   Thermocouple    in   Non- Isothermal   Flow 
The   behavior    of    a   thermocouple    over    a    fairly    large   range    of 
temperature    may   be    expressed    as 


E   =    E      +   aT    . 
o 

The  following  fluctuating  quantities  are  introduced: 


(19) 


E   =    E   +    e'  and  T    =    T   +    t'     . 

Substitution    of  Equations    20    into    19   yields 


(20) 


E+e'=E     +  a(T  +  t') 
o 


(21) 


Time  averaging  of  Equation  21  yields 

E  =  E   +  a?   .  (22) 

o 


Subtracting  Equation  22  from  21  yields 

(23) 

Squaring,  time  averaging,  and  taking  the  square  root  of  Equation  23 

gives  — 

/  e'2   =  a/  t'2   .  (24) 

Hot-Wire  in  Non-Isothermal  Flow 

4 
The  following  parallels   Corrsin's   treatment  lor  a  constant 

current  hot  wire  anemometer. 

Equation  8  is  still  valid,  however,  P  and  k  are  no  longer  constants 

Substitution  of  Equation  11  into  8  yields 

E2  _    (lL_  Mk    ,  rTr£D1/2      r    k   ,1/2        1/2    1/2,  1/2.    1/2 

R  (R  -R  )"  (aR  A)k+(   aR    '[Cm     B)p    u    k    Cp    '  (2D 

w  w  a      o  op  r 


It  is  assumed  that 

k 


C  u   Pr 
P 

is  constant,  and  that  C   itself  is  also  constant. 

P 

The  following  mean  and  deviation  quantities  are  considered: 


Ra  =  Ra  +  V  P=P  +  P' 


(26) 


k=k+k  E=E+e' 

u   =  u  +  u' 


Substitution  of  Equations  26  into  2!  yields 

(i+f')2   ■  ffAKk*)  + 

R  (R  -R  -r  ')    U  o  . 

wwaa  .i/l 

+  (^i—  .t^]   B)(-+,')1/2(u+u')1/2(k+k')1/2 
o 


(27) 


According  to  the  Binomial  Theorem, 


(E  +    e')2 

1 

R    (R      -    R    )     ' 
ww          a 

1- 

r    ' 

a 

R   -  R 
w      a 

(E+e')2 

R  (R  -R  ) 

w  w  a 


(1  +  —2—  ) 


R  -R 

w  a 


(28) 


Note  that 


(E  =  e')2  -  E2  +  2Ee'  . 


(29) 


Using  Equations  28  and  29  for  the  left  side  of  27,  and  expanding 
the  right  side  of  27,  and  neglecting  small  quantities  yields: 


E              a         2  E  e '      ,   tt  'I    ,  .  .—   ,  , . 
__  +  __  +  __  =  (      A)(k  +k.)  + 

R  (R  -R  )    R  (R  -R  )     R  (R  -R  )        o 


w  w   a     w  w  a 


w  w  a 


(30) 


1/2  C   1/2   _ 

+  (-1— h^"]    B)(Puk  +  Duk'  +  pu'k  +  {.   'u  k) 

tR  Pr 
o 


1/2 


But 


(puk  +  puk1  +  pu'k  +  p *uk) 


1/2 


1/2,,  ,  lrk' 


2  k    u    p 


(31) 


Thus  Equation  30  becomes 

—2 
—2  E   r    '  — 

L"3-  +  ^—J+        2Ee'-        "    <  -JT  A)(k  +   k')    +  (32) 

R    (R   -R    )       R    (R   -R    )  R    (R   -R   )  o 

wwa  wwa  wwa 


»         ..  >         „  • 


„_l/2      C      1/2  . 

+  (^f—    [^]         B)(p^k)1/2(l+±    [f-  +  ^-  +  -         ]). 

aR  Pr  2      k  u  p 

o 

The  time  average  form  of  Equation  32  is 

72  o  on1/2    C    1/2        1/9 

v    rv     -  =  (^~A)k+  (  ~^~  [  77  ]    B)(puk)    •  (33) 

R  (R  -R  )      o  o 

wwa 


Subtraction  of  Equation  33  from  32  gives 

I2    r    ' 

a  ,  2Ee'  ,    rr  *  , 

=— 2   +  =~  =    (  ^R~  A)k      +  (34) 

R   (R  -R   )  R    (R  -R   )  o 

wwa  wwa 


0_l/2      C      1/2        — -1/2      .  , 
o 


The  following  relationships  are  known: 


'      —  t' 
k'  =  nt'       p   =  -  p 


T  (35) 

a 

Thus  Equation  34  becomes 

E2  R   a  - 

o  .  2E  .  .    nTifc      .      , 

___   ti    +  __  e'    =    (  — ~  A)t*    + 

R    (R  -R   )  R    (R  -R   )  o  H6) 

wwa  wwa  \JU' 

onl/2   C   1/2    — -p.  1/2 
rwD      p       (puk)     .   t1    t1    u'  . 

+  (  aR     [Tr]         B)     2 (n  V  ~  T  +  ~  }  ' 

o  a    u 


i-   '     t-  '  ' 

Since   n  —  -  —  <<  —  >  Equation  36  becomes 

k    T        u 
a 


10 


2E 


R    (R   -R    ) 
•  ww     a 


,     .  Wnl/2c    1/2L,1/2k1/2  " 


e      = 


.   2a,      Pr1/2   5    1/2 
o 


u'    + 


[■ 


aR 


-2 

E      R     a 


o  R    (R    -R 

ww     a 
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Equation  37  may  be  written 

rR  (R-R  ),€D1/2C  1/2BP1/2R1/2 
w  w  a  p 

~~-  1/2  -  1/2 

4  E  a  R   Pr    u 
o 


u'  + 


+ 


R  (R  -R  )  nir£A 
w  w   a 

2  E  a  R 


E  R  a 
o 

2  (R  -R  )  J 
w   a 


t* 


(38) 


The  coefficient  of  u'  represents  the  sensitivity  to  velocity  fluctuations, 
and  the  coefficient  of  t'  represents  the  sensitivity  to  temperature 
fluctuations.   These  sensitivities  could  in  principle  be  calculated, 
since  all  factors  are  either  measurable  quantities,  physical  properties, 
or  known  constants.   In  practice,  however,  this  calculation  yields  values 
of  the  sensitivities  that  do  not  agree  with  those  determined  by  actual 
calibration.   This  discrepancy  is  due  to  two  of  the  supposedly  known 
quantities  in  Equation  38,  namely  <y   and  B.   King's  value  for  B,  since 
it  is  based  on  his  potential  flow  assumption,  is  seriously  in  error.   The 
value  of  ot      found  in  tables  of  physical  properties  may  also  be  in  error 
because  the  value  of  a      applicable  to  the  wire  in  use  depends  very  strongly 
on  the  history  of  that  wire.   In  addition,  contamination  of  the  wire 
changes  the  value  of  D,  the  wire  diameter,  and  the  wire  length  I   is  some 
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effective  length  rather  than  the  actual  length.   Thus  it  is  impossible 
to  accurately  determine  the  sensitivities  by  direct  calculation.   In- 
stead the  sensitivities  are  determined  by  calibration. 

Although  the  expressions  for  the  sensitivities  in  Equation  38  do 
not  provide  quantitative  information,  they  do  indicate  the  direction  of 
variation  in  the  sensitivities  as  various  parameters  are  changed.   Thus, 
the  velocity  sensitivity  should  decrease  with  increasing  velocity,  but 
velocity  should  not  affect  the  temperature  sensitivity  to  first  order. 
The  velocity  sensitivity  should  increase  with  increasing  wire  temperature 
(resistance) . 

In  order  to  determine  the  effect  of  parameter  variation  on  the 
temperature  sensitivity,  further  consideration  of  the  two  terms  that 
comprise  it  is  necessary.   Equation  8  shows  that,  all  other  parameters 
remaining  constant,  an  increase  in  T   (i.e.  a  positive  t')  requires  a 

9. 

2 
decrease  in  E   (i.e.  a  negative  e')  in  order  to  maintain  the  equality. 

Thus,  the  overall  sign  of  the  temperature  sensitivity  must  be  negative. 

Since  all  factors  are  inherently  positive,  this  implies  that 


E   R    a                  R    (R    -R    )    niri  A 
2 >         w      w      a ..  (39) 


2(R    -R    )  2   E  QR 

w      a  o 

ER    a 

™,  .  .  o  . 

Thus    an    increase    in  tttt r— .  causes    an    increase    in 

t-  ( R  -  R  ) 
w   o 

the  magnitude  of  the  temperature  sensitivity,  and  an  increase  in 

R    (R    -R    )nTT0A 
w     w      a         ^ 

tt= — causes    a   decrease. 

2E  a  R  S2 

o  E 


Equation  33  shows  that  the  quantity     R  (R  -R  ) 

w  w   a 

is  independent  of  T   -  f  .   That  is  if  the  parameter  T   -  T   is  varied 

w    a  w    a 
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-2  -  -2 
E           remains  constant.   Further  if  R   is  constant,  E   is 
a 

R  (R   -  R  ) 
ww    a 

2 
proportional    to  R    .        Now 

w 


E    R     o  R  -2  ^   R 

o  =     _w  E  o       , 


(40) 


2    (R    -    R    )         E  R    (R    -R    )       2 

w        a  wwa 


and 


R    (R    -R    )    nirtih         -      R    (R    -R    )       n*/?A 
wwa  X/=E        ww      a  I 


2   E    ot  R 
o 

E2 

R   a 
o 

!    R 

w 

-2 
and        E 

are    in< 

E 

R    (R    -R    ) 

wwa 

(41) 


are    independent    of  T    -T    ,    it 

w  a 


follows  from  Equation  40  that  variation  of  T  -T   has  no  effect  on 

w      a 

E   R    a 
o 


2(o      d    \  to    first    order.      Similarly    it    follows    from 

V  w"    aJ  R    (R   -R    )    n*/A 

ww      a  " 

Equation   41    that  —  increases    linearly  with 

2  Ea  R 
o 

increasing  T  -T  .   As  previously  shown,  an  increase  in  this  term  causes 
w  a 

a  decrease  in  temperature  sensitivity.   Thus  an  increase  in  T  -T   causes 

w      a 

a  decrease    in    temperature    sensitivity.      Thus    it    appears    that   velocity 
fluctuations    could   be    measured   directly   by    operating   at   ver)    high 
temperature   differences    and    that    temperature    fluctuations    could   be 
measured   directly   by    operating   at    very    low   temperature   differences. 
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If  consideration  is  given  to  changing  the  dimensions  and  material 
of  the  wire,  further  effects  on  the  sensitivities  occur.   Velocity 
sensitivity  increases  with  wire  length  while  the  temperature  sensitivity 
decreases.   An  increase  in  wire  diameter  causes  an  increase  in  velocity 
sensitivity  and  no  change  in  temperature  sensitivity  to  first  order.   An 
increase  in  temperature  coefficient  causes  a  decrease  in  velocity 
sensitivity  and  an  increase  in  temperature  sensitivity. 

Once  the  desired  information  on  trends  in  sensitivity  is  recognized, 
it  is  convenient  to  resort  to  a  more  compact  notation: 

R    (R   -R    )^D1/2C1/2B    V.ll2im 

ww      a       ** p_ 

Su   5  "      4  EC,  R      P    1/2   C,1/2 
o      r 

(42) 


and 


ERff  R  (    R  -R    )n   /An 

q      =         o  -      _w      w      a        *- 

2(R    -R    )  2EaR 

w      a  o 


(43) 


thus  Equation  38  becomes 


e'  =  Su  u'  -  ST  t'.  (44) 


Squaring  and  time  averaging  yields 


2     2  2    2 

e,Z  =  S    -  2S   S   u't'  +  S    t'  .  (45) 

u       u   T         T 

The  sensitivities  S   and  Sm  must  be  determined  from  calibration 

u      T  

2 
This  is  discussed  in  the  section  on  procedure.   The  quantity,  e1  , 
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2 
is  experimentally  measured.   From  known  values  of  e1  ,  S  ,  and  S  ,  a  set 

of  simultaneous  equations  of  the  form  of  Equation  39  can  be  formed. 


This  set  can  then  be  solved  for  u'  ,  u't1,  and  t'  .   This  solution  is 
discussed  in  the  RESULTS  section. 
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APPARATUS 

Flow  System 

The  flow  system  is  shown  schematically  in  Figure  1.   A  Craftsman 
Model  315.16970,  3HP  industrial  vacuum  cleaner  was  operated  as  a  blower 
to  supply  air  to  the  system.   A  three  inch  galvanized  duct  carried  the 
air  to  the  highest  point  of  the  flow  system  where  the  duct  tapered  to 
meet  the  1.0625  inch  inside  diameter  Type  M  copper  tubing  that  comprised 
the  remainder  of  the  loop.   The  three- inch  duct  was  fitted  with  a  bleed 
line  and  damper.   The  flow  rate  in  the  loop  could  be  varied  by  using  the 
damper  to  alter  the  division  of  air  between  loop  and  bleed  line.   Above 
the  bleed  line  in  the  duct  was  an  orifice  meter  for  measuring  air  flow 
rate  in  the  loop.   (The  calibration  curve  for  the  orifice  meter  is 
given  in  Appendix  B.)   In  the  copper  tube,  the  air  passed  downward 
through  an  83  inch  unheated  velocity  developing  section.   This  was 
followed  by  a  58  inch  heating  section  which  was  wound  with  nichrome 
heating  wire  and  covered  with  1.125  inches  of  pre-formed  pipe  insulation. 
Probes  were  inserted  upward  into  the  copper  tube  as  shown.   All  measure- 
ments were  made  in  the  end  plane  of  the  heating  section. 

Iron-constantan  bulk  thermocouples  were  installed  at  the  beginning 
of  the  velocity  developing  section  and  at  the  end  of  the  heating  section, 
and  wall  thermocouples  were  installed  at  the  end  of  the  heating  section. 
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Figure  1.   Flow  Diagram  of  Apparatus 
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Probes  and  Traversing  Mechanisms 

The  rapid  response  thermocouple  used  (a  Heat  Technology  Laboratory 
TCFW  202-ChA-5)  is  shown  in  Figure  2.   The  response  time  of  this  thermo- 
couple has  been  estimated  by  Rodriguez-Ramirez   to  be  5  milli- seconds . 

The  hot-wire  probe  used  (a  Thermo  Systems  1210-T1.5)  is  shown  in 
Figure  3. 

Positioning  of  the  rapid  response  thermocouple  was  by  means  of  the 
traversing  mechanism  shown  in  Figure  4.   When  in  place  as  shown  in  Figure 
1  the  mechanism  was,  of  course,  inverted.   The  ball  joint  permitted 
radial  movement  of  the  probe  across  the  tube.   Radial  position  was 
measured  with  the  dial  guage  shown. 

Figure  5  shows  the  traversing  mechanism  used  to  position  the  hot- 
wire in  the  loop.   This  traversing  mechanism  replaced  the  thermocouple 
traversing  mechanism  at  the  bottom  of  the  loop  when  hot-wire  measure- 
ments were  being  made.   The  upper  block  in  Figure  5  could  be  moved 
laterally  so  as  to  place  the  hot-wire  probe  at  the  desired  radial  position, 
while  the  lower  plate  shown  in  Figure  5  was  bolted  to  the  loop.   A  pitot 
tube  and  a  local  mean  thermocouple  were  inserted  through  the  auxiliary 
probe  entrances. 
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Figure  2.   Fast  Response  Thermocouple 
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Figure    5.      Hot-Wire  Traversing  Mechanism 
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Inst rumen tat  ion 

Instrumentation  for  hot-wire  measurements  consisted  of  a  Thermo 
Systems  Type  1010  anemometer,  and  a  Ballantine  Model  320A  True  Root- 
Mean-Square  Electronic  Voltmeter. 

For  rapid-response  thermocouple  measurements,  instrumentation  consist- 
ed of  a  Tektronix  122  Amplifier  (See  Appendix  C  for  frequency  response.) 
and  the  Ballantine  Voltmeter. 

Bulk  and  wall  thermocouple  measurements  were  made  with  a  Dana  Model 
2850  D.C.  amplifier  and  a  Beckman  Model  4011R  digital  voltmeter. 
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PROCEDURE 


Hot-wire  in  Isothermal  Flow 


The  first  step  was  calibration.   The  constant  C  defined  by  Equation 

17  and  appearing  in  Equation  18  can  be  determined  from  mean  voltage  and 

-2 
velocity  measurements.   It  is  seen  from  Equation  14  that  if  E 


R  (R  -R  ) 
w  w  a 


were  plotted  as  a  function  of  u~,  the  result  should  be  a  straight  line 
with  slope  C. 

The  Calibration  of  the  hot-wire  probe  was  carried  out  with  the  probe 
in  place  in  the  loop  and  located  at  the  tube  centerline.   Variation  of 
u  was  achieved  by  varying  the  flow  rate  in  the  loop  using  the  bleed  line 
damper.   The  quantity  E  was  measured  at  seven  or  eight  different  flow 
rates.   The  value  of  u  at  the  centerline  corresponding  to  each  flow  rate 
was  measured  with  a  pitot  tube.   The  data  were  plotted  as  indicated  above, 
a  straight  line  was  fitted  to  the  points,  and  the  slope  was  measured. 
Calibration  curves  are  shown  in  Appendix  E. 

The  second  step  was  the  making  of  fluctuation  measurements.   Equation 
18  gives  the  relationship  between  the  rms  velocity  fluctuation  and  the 
rms  voltage  fluctuation. 

The  flow  system  was  operated  at  Reynolds  Numbers  of  30,700,  41,500, 
and  49,800.   For  each  Reynolds  Number,  measurements  were  made 


ofVe'2 


at    eight    points   between    the    center   and    the  wall    of    the    tube    corresponding 
to   r/R    of    0,    0.125,    0.250,    0.375,    0.500,    0.625,    0.750   and    0.875.      A 
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pitot  tube  was  used  to  measure  u  at  the  same  points. 

After  all  fluctuation  measurements  had  been  made,  the  calibration 
procedure  was  repeated  to  make  certain  that  the  constant  C  had  not 
changed  during  the  run. 

Rapid  Response  Thermocouple  in  Non-Isothermal  Flow 

The  flow  was  adjusted  so  that  the  Reynolds  Numbers  in  this  part  were 

the  same  as  those  in  the  previous  section.   The  applied  heating  was  430 

watts  in  all  cases. 


Measurement  was  made  of  \  e '   from  the  rapid  response  thermocouple 
as  a  function  of  radial  position  for  each  Reynolds  Number.   No  calibration 

was  done,  since  the  constant,  a,  in  Equation  24  was  assumed  to  be  that 

13 
found  in  standard  thermocouple  tables   .   This  assumption  is  of  question- 
able validity. 

A  conventional  thermocouple  was  used  to  measure  the  centerline 
temperature,  and  the  wall  thermocouples  were  used  to  measure  the  wall 
temperature  for  each  Reynolds  Number. 

Hot-Wire  in  Non- Isothermal  Flow 

The  procedure  of  calibration  and  measurement  outlined  below  repre- 
sents one  of  the  contributions  of  this  work.   It  is  anticipated  that 
this  procedure  will  be  used  for  the  glycol  and  mercury  measurements 
previously  mentioned. 

As  shown  in  the  THEORETICAL  BACKGROUND  section,  the  sensitivities 

S   and  S   depend  on  two  parameters:   u  and  T  -T  .   Thus,  in  calibration 
u      T  w   a 

it  is  necessary  to  determine  the  functional  dependence  on  each  parameter. 
Calibration  for  temperature  sensitivity  was  also  carried  out  with 
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the  probe  in  place  and  at  the  tube  center  line.   Velocity  variation  was 
again  achieved  by  use  of  the  bleeder  line  damper,  and  u  was  measured  with 
a  pitot  tube.   In  this  way  five  different  values  of  u  were  used,  and 

for  each  value  of  u  thirteen  values  of  T  -T  were  used. 

w   a 

These  data  were  plotted  as  E  versus  u  for  constant  values  of  T  -f 

w    a 

and  as  E  versus  T  -T   for  constant  values  of  u.   These  curves  are  shown 
w    a 

in  Appendix  G.   The  sensitivities  S   and  Sm  are  the  derivatives  of  these 

rr  U         T 

curves.   These  derivatives  were  obtained  from  analytical  differentiation 

of  parabolas  fitted  to  the  data  by  the  method  of  least  squares.   The 

selection  of  parabolas  is  arbitrary,  but  the  relatively  slow  and  smooth 

variation  of  the  data  makes  a  parabola  a  satisfactory  approximation. 

Note  that  T  -T   must  be  computed  from  measurements  made  of  R  -R  . 
w   a  w   a 

The  resistance- temperature  characteristic  of  the  hot-wire  probe  was 

experimentally  determined  to  be  as  shown  in  Appendix  D. 

Once  S   and  S   were  known  as  functions  of  u  and  T  -  T  ,  fluctuation 
u      T  w    a 

measurements  were  made.   The  rms  voltage  fluctuation  was  measured  at 
r/R  =  0  and  r/R  =  .5  for  each  flow  rate  and  temperature  difference  con- 
sidered in  calibration. 

When  all  fluctuation  measurements  had  been  completed,  the  calibration 
was  repeated  to  determine  the  drift  in  S   and  S  . 
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RESULTS 


Hot-Wire    in    Isothermal   Flew 

Figure    6    shows    the    axial    turbulent    intensity   at    the    centerline    of 
the    pipe    compared   with    the   work    of   Sandborn. 

Figure    7    shows    axial    turbulent    intensity   as    a    function   of    radial 
position    for    three   different   Reynolds   Numbers    compared  with    the   work   of 

T  f  10 

Laufer . 

The  agreement  of  present  results  with  those  of  Sandborn    and 
Laufer   is  rather  good.   The  data  used  by  Sandborn  to  obtain  his 
correlation  show  a  scatter  of  j  10  percent.   The  results  of  the  present 
work  are  well  within  this  range. 

Rapid  Response  Thermocouple  in 
Non- Isothermal  Flow 

Figure  8  shows  the  temperature  fluctuation  intensity  as  a  function 

14  15 

of  radial  position.   Rodriguez-Ramirez   ,  and  also  Tanimoto  and  Hanratty 

have  published  results  of  this  nature.   The  present  results  are  for  a 
range  of  Reynolds  Numbers  higher  than  the  range  considered  by  Rodriguez- 
Ramirez;  however,  a  preliminary  run  was  made  in  which  the  flow  rate  and 
heat  input  were  the  same  as  those  in  one  of  Rodriguez-Ramirez' 
runs.   The  results  of  this  run  deviated  from  those  of  Rodriguez-Ramirez 
by  less  than  5  percent.   The  temperature  gradients  used  by  Tanimoto  and 
Hanratty  were  so  much  larger  than  those  used  in  the  present  work  that 
comparison  is  not  possible. 
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Figure  8.   Temperature  Turbulence  Intensity  in  Air 
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Hot-Wire  in  Non- Isothermal  Flow 

Figures  9  and  10  show  the  dependence  of  the  sensitivities  S   and  S 

b  r  u      T 

on  the  parameters  u  and  T  -T  .   The  velocity  sensitivity,  S  ,  increases 

w   a  u 

with  increasing  temperature  difference  and  decreases  with  increasing 
velocity  as  predicted  from  Equation  38.   The  temperature  sensitivity 
decreases  with  increasing  temperature  difference  as  predicted  from  Equation 
38.   Note  that  there  is  a  secondary  dependence  on  velocity  not  predicted 
by  Equation  38.   This  dependence  is  due  to  the  velocity  dependence  of  the 
effective  wire  length. 

The  investigation  of  various  calculating  schemes  including  the 
selection  of  the  best  one  for  analyzing  the  non- isothermal  data  was  an 
important  part  of  the  present  work.   As  pointed  out  in  the  section  on 
theoretical  development,  Equation  45  gives  the  relationship  between  the 


unknown  quantities  u'  ,  u't'  and  t'   and  the  known  quantities  S  ,  S  , 

2 
and  e1   .   Since  there  are  three  unknowns,  at  least  three  equations  of  the 

form  of  Equation  39  are  needed.   These  can  be  obtained  by  operating  the 

hot-wire  at  least  three  values  of  T   (and  hence  of  T  -T  ) .   For  each 

w  w   a 

value  of  T  -T  ,  the  values  of  S   and  S   appropriate  to  the  mean  velocity, 
w   a  u      T  

2 
u,  of  the  flow  are  read  from  Figures  9  and  10,  and  e1   is  measured.   Thus, 

for  each  value  of  T  -T   an  equation  is  formed. 

w   a 

Once  the  set  of  equations  is  obtained,  a  number  of  methods  of  solution 
are  available.   The  most  obvious  one  is  the  simultaneous  solution  of 
three  such  equations  by  a  method  such  as  Gaussian  reduction.     However, 
the  numerical  calculations  involve  the  subtraction  of  large  numbers  whose 
difference  is  small,  which  has  the  effect  of  greatly  magnifying  the  errors 

in  the  data.   Since  the  errors  in  S   and  S   were  already  large,  the  result 

2  2 

was  errors  in  the  calculated  values  of  u1  ,  u't',  and  t1   of  the 
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Figure    9.      Hot-Wire   Temperature   Sensitivity 
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Figure    10.      Hot-Wire   Velocity   Sensitivity 
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order  of  500  percent . 

The  Kovasznay  Fluctuation  Diagram    is  a  graphical  method  of  solving 

2 
the  set  of  equations  in  which  equation  39  is  divided  by  S    to  give: 


^4=  (  -^  )2  u'2  -  2(  ^  )  u't'  +  t'2  .  (46) 

S  J  bT  ^T 

T 


P. 2  S 

e  u 

Then  the  quantity    J  is  plotted  versus  ■—  .   Note  that  analytically 

S^  T 

T 
the  relationship  between  the  two  is  described  by  a  parabola.   By  the 

method  of  least  squares  a  parabola  is  fitted  to  the  data.   The  computer 

program  used  to  do  this  fitting  is  listed  in  Appendix  F.   The  intercept 

"T2  — t  S 

e  2  u 

of  this  parabola  with  the axis  is  t1    ,  the  slope  at   —  =0is 

ST  ^  T 

-  2u't'  ,  and  the  curvature  of  the  parabola  is   u'   .   This  method  was 

S 

used  on  the  present  data,  but  the  need  to  extrapolate  the  data  from  — -  of 

S  S     ST 

u  u 

some  non-zero  value  to   Q   =0  and  then  determine  the  slope  at  — —   =  0 

T         T 

resulted  in  a  lar^e  uncertainty  in  u't'  .   This  same  problem  was  encountered 

9 

by  Arya  and  Plate  ,  who  suggested  that  if  an  independent  measurement  of 

~2 

t     could  be  obtained  to  locate  the  intercept  of  the  curve  accurately, 

2 

then  reliable  values  of  the  remaining  two  unknowns   u't'  and   u'   could 

be  determined  by  the  Kovasznay  method.   This  procedure  was  followed  with 

the  present  data.   The  value  of   t    used  was  the  value  found  using  the 

rapid  response  thermocouple.   The  results  were  still  rather  uncertain, 

but  the  uncertainty  was  less  than  for  any  of  the  other  methods.   (The 

Kovasznay  Fluctuation  Diagrams  for  representative  cases  arc  shown  in 

Appendix  H.)  

2  — 
The  numerical  values  of  *    u'  /u  c   presented  in  Figures  11  and  12 

deviate  seriously  from  the  isothermal  results  wilIi  which  it  is  thought 

they  should  agree.   (See,  for  instance  Corrsin  and  Uberoi   )  The  trends 
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of  the  results  with  Reynolds  Number  and  radial  position  are  in  accordance 
with  expectations. 


The  calculated  values  of  u't'/u.  (T  -T  )  are  presented  in  Figure  13 

tot  ° 

Turbulence  at  the  centerline  in  pipe  flow  for  the  Reynolds  numbers  used 


in  the  present  work  should  be  isotropic.   Thus  u't'  and  u't'/u   (T  -T  ) 

tot, 

should  both  be  zero  at  the  centerline.   The  centerline  results  in 
Figure  13  do  not  exhibit  this  behavior.   However  the  error  in  these  re- 
sults is  roughly  estimated  in  Appendix  I  as  1007o. 

Away  from  the  centerline  the  turbulence  should  no  longer  be  iso- 


tropic, and  u't'  should  no  longer  be  zero.   With  heating  from  the  wall 


as  in  the  present  case,  u't*,  and  hence  u't'/u   (T  -T  )  should  be  negative 

cot 


This  trend  is  followed  by  the  present  results  in  that  u't'/u  (T  -T  )  at 

tot 

r/R  =  .5  has  a  larger  magnitude  than  at  r/R  =  0  and  in  that  it  has  a 


negative  sign.   The  quantity  u't'/u  (T  -T  )  cannot  continue  to  increase 

tot 

indefinitely  as  Reynolds  Number  increases,  for  this  would  mean  no  dis- 
sipation was  taking  place.   The  increasing  trend  in  the  results  shown  in 
Figure  13  for  r/R  =  .5  is  therefore  not  anticipated,  but  it  is  possible 

so  long  as  the  trend  reverses  itself  at  higher  Reynolds  Numbers.   Arya 

9 
and  Plate   found  that  in  the  boundary  layer  growing  on  the  wall  of  a  wind 


tunnel  u't'/u  (T  -T  )remained  nearly  constant  with  increasing  Reynolds 
Number.   Again  with  an  estimated  error  of  100  percent   little  significance 
can  be  attached  to  present  results. 
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DISCUSSION 


The   numerical   results    for\    u'    /u      and   u't'/uT 

t  L 

are  in  general,  too  large.   In  view  of  the  way  in  which  these  intensities 

are  obtained  from  the  Kovasznay  Fluctuation  Diagram,  this  suggests  that 

S   is  too  small.   The  temperature  sensitivity  S_  affects  both  coordinates 
u  T 

in  the  same  way,  and  therefore,  errors  in  S  would  be  self-compensating. 

S 

However,  errors  in  S   affect  only  the  —  coordinate.   Too  small  a  value 

u  Sm 

T 

of  S  would  cause  too  large  a  slope  and  curvature  in  the  fitted  parabola 
u 


and  too  large  a  value  of  u'   and  u't'. 

Because  of  the  long  time  required  for  the  system  to  come  to  thermal 

equilibrium  after  each  change  in  flow  rate,  the  data  used  to  compute  S 

for  a  given  value  of  T  -T   were  taken  over  a  period  of  approximately 

w   a 

eight  hours.   It  is  reasonable  to  expect  an  increase  in  the  thickness  of 
the  layer  of  contamination  on  the  wire  during  this  period  since  the  wire 

was  continuously  exposed  to  an  inadequately  filtered  flow. 

3 
Kronauer   has  shown  that  dirt  accumulations  on  fine  wires  at  low 

flow  rates  increase  the  heat  transfer  rate  because  the  heat  transfer  area 

increases  more  rapidly  than  the  heat  transfer  resistance  as  dirt 

accumulates.   Thus,  a  dirty  wire  would  require  a  higher  value  of  E  than 

a  clean  one  to  maintain  it  at  a  given  temperature. 
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The  wire  used  for  the  present  measurements  was  examined  under  a 
microscope,  and  a  significant  amount  of  dirt  was  observed  on  it. 

In  the  present  calibration,  E  was  measured  at  the  highest  calibration 
velocity  first,  and  the  velocity  was  decreased  step  by  step  until  the 
lowest  value  was  reached.   As  anticipated,  E  decreased  with  decreasing 
velocity,  but  the  increasing  layer  of  contamination  apparently  worked 
counter  to  this  trend.   Thus,  in  the  calibration  procedure  used  the  rate 
of  decrease  of  E  was  slower  than  if  contamination  had  not  been  a  problem. 
The  result  was  a  value  of  S   that  was  too  small. 

su 

Larger  values  of  S   would  result  in  larger  values  of  —  at  the  same 

m 

2     2 

value  of  e '   -  t    on  the  Kovasznay  Fluctuation  Diagram  and  thus,  a 

smaller  slope  and  curvature  of  the  fitted  parabola.   The  calculated  values 


of  \  u '    and  u't'  would  then  be  smaller  and  in  better  agreement  with  the 
isothermal  results. 


It  was  not  possible  to  determine  the  errors  in  \  u'   and  u't' 
caused  by  the  drift  in  E  noted  above  because  it  was  not  possible  to 
measure  what  E  would  have  been  without  contamination.   Thus  it  was 
necessary  to  use  a  different  method  to  estimate  these  errors.   The  method 
used  is  outlined  in  Appendix  I,  and  the  error  estimates  are  36  percent 


in  \  u  '    and  100  percent  in  u't'.  As  noted  in  Appendix  I,  these  figures 
do  not  represent  true  confidence  limits  on  the  results.   They  are  only 
rough  estimates  of  the  errors  present. 

The  calibration  curves  for  before  and  after  as  shown  in  Figures  24 
through  26  show  such  a  large  drift  that  all  data  taken  in  between  should 
be  discarded.   This  was  not  done,  since  a  number  of  attempts  to  obtain 
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better  results  failed.   Note  then,  that  all  fluctuation  intensities  ob- 
tained from  the  present  non- isothermal  hot-wire  measurements  are  of  very 
doubtful  reliability. 

Due  to  inexperience  in  this  area,  a  number  of  deviations  from  good 
experimental  practice  were  present.   First,  although  the  wire  was 
examined  under  a  microscope  after  all  measurements  had  been  made,  no 
examination  was  made  before  or  during  experimental  runs.   Thus  the 
hypothesis  of  an  increasing  dirt  layer  cannot  be  verified.   Second,  no 
attempt  was  made  to  clean  the  wire  during  the  course  of  experimental 
runs.   Because  of  these  errors  it  is  not  possible  to  analyze  in  detail 
or  to  correct  for  the  drift  that  occurred. 

In  future  work,  provision  should  be  made  for  careful  filtration  of 
the  fluid  and  a  means  of  rapid  calibration  to  minimize  drift.   Two  cali- 
brations, one  with  increasing  velocity  and  one  with  decreasing  velocity, 
should  be  made  to  check  for  contamination.   If  it  is  still  not  possible 
to  reduce  drift  to  an  acceptable  level,  provision  should  be  made  for  the 
complete  withdrawal  of  the  probe  from  the  flow. 

A  problem  encountered  with  the  present  apparatus  was  difficulty  in 
maintaining  flow  rate  and  heat  input  at  desired  levels.   Long  term 
temperature  variations  of  two  or  three  degrees  persisted  even  after  the 
system  had  come  to  "equilibrium".   This  was  apparently  caused  by  cycling 
of  the  room  temperature. 

The  nature  of  the  problems  encountered  in  the  present  work  serve  as 
a  guide  in  planning  additional  steps  in  the  overall  research  program. 
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SUMMARY  AND  CONCLUSIONS 

This  work  has  developed  procedures  for  the  calibration  of  a  hot- 
wire in  non-isothermal  flow.   Further,  the  various  calculation  schemes 

available  for  analyzing  non- isothermal  flow  fluctuation  data  have  been 

9 
evaluated  and  the  Arya  and  Plate   modification  of  the  Kovasznay  Fluctuation 

Diagram   method  has  been  found  to  be  the  most  satisfactory. 

Apparatus  limitations  and  deviations  from  good  experimental  practice 
precluded  the  production  of  reliable  numerical  results. 

The  findings  of  this  investigation  are  an  important  preliminary  step 
in  a  comprehensive  research  program  of  non- isothermal  flow  measurements 
with  hot-wire  and  hot- film  anemometers. 
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APPENDIX  A 

TABLES    OF    EXPERIMENTAL  DATA 

Table    1 

Orifice   Meter   Calibration 

AP   Orifice  ;.  P    Pilot    Tube 

Meter      (in)*  At  (in)* 

4.553  2.173 

3.151  1.595 

2.172  1.123 

1.114  0.571 

*   Man  ome  t  e  r    Fluid  S  PC  .    0.797 
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Table  2 
Tektronix  Amplifier  Frequency  Response 


Frequency 

Vc 

1 1  a ge  i  n 

Voltage  Out 

(C  PS  ) 

(Mi 

llivolts) 

(Volts) 

5 

.200 

.200 

7 

.200 

.214 

9 

.200 

.222 

20 

.200 

.238 

40 

.200 

.244 

60 

.200 

.245 

80 

.200 

.245 

100 

.200 

.244 

200 

.200 

.239 

300 

.200 

2.31 

400 

.200 

.226 

500 

.200 

.217 

600 

.200 

.210 

700 

.200 

.200 

800 

.200 

.193 

900 

.200 

.185 

1000 

.200 

.  177 

1100 

.200 

.170 
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Table    3 
Calibration    for    Isothermal  Re=49,800  Run 


AP  (  in   of  — - 

Pitot   Tube   SPG.  797        R    (OHM)  R    (OHM)  E  (VOLTS )              /  e'    (MILLIVOLTS) 
Fluid) 

3.890  7.57  10.60  3.07  13.0 

3.394  7.54  10.56  3.04  13.0 

2.948  7.52  10.53  2.98  12.9 

2.407  7.50  10.50  2.93  12.8 

1.935  7.48  10.47  2.88  12.7 

1.369  7.46  10.44  2.81  12.5 

0.844  7.43  10.40  2.71  12.5 

0.330  7.38  10.33  2.54  12.6 

Before     t 

4.004  7.57  10.60  3.08  

0.737  7.42  10.39  2.70  

After       t 
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Table  4 
Data  for  Isothermal  Re=49,800  Run 


Radial  /\P                      (in    of                                                 - 

Position  (in)                         Pitot   Tube      SPG. 797  L(Volts)  /  e'    (Millivolts) 
Fluid) 

0.0000  3.984  3.07  12.9 

0.0625  3.964  3.07  12.9 

0.1250  3.841  3.07  16.5 

0.1875  3.691  3.06  19.5 

0.2500  3.433  3.04  23.6 

0.3125  3.184  3.02  27.5 

0.3750  2.870  2.97  31.7 

0.4375  2.411  2.94  36.4 


Table    5 
Calibration    for    isothermal   Re=30,700   and   41,500  Runs 


AP  Pitot   Tube 

(in.    of    SPC.797   Fluid) 


R  (Ohms) 
a. 

R  (Ohms) 
w 

E  (Volts) 

2 
v    e   (Millivolts) 

7.55 

10.57 

3.05 

12.8 

7  .V} 

10.54 

2.98 

12.6 

7.51 

10.51 

2.94 

12.5 

7.48 

10.47 

2.88 

12.4 

7.46 

10.44 

2.82 

12.4 

7.44 

10.42 

2.75 

12.3 

7.40 

10.36 

2.61 

12.4 

7.52 

10.53 

2.96 

—  —  —  — 

7.44 

10.42 

2.75 

3.533 
2.999 
2.500 
2.000 
1.496 
1.020 
0.500 

Before  f 

2.780 
1.020 

After     t 
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Table  6 
Data  for  Isothermal  Re=30,700  Runs 


Radial         ,.       . 
(in .  ) 
Position 

AP 
(in 

Pitot   Tube 
,.  SPG.  797   FT 

aid) 

E 

(Volts) 

/   e'2 (Millivolts) 

0.0000 

1.480 

2.81 

12.2 

0.0625 

1.460 

2.81 

12.9 

0.1250 

1.400 

2.80 

15.1 

0.1875 

1.330 

2.79 

18.2 

0.2500 

1.240 

2.78 

21.9 

0.3125 

1.130 

2.76 

24.8 

0.3750 

1.02  0 

2.74 

28.6 

0.4375 

0.850 

2.71 

32.5 

Table  7 
Data  for  Isothermal  Re=4l,500  Run 


Radial          ..       . 

■a    -•<-•             (in-  > 

Position 

AP 
(in 

Pitot   Tu 
.SPG797 

be 
Flui 

d) 

E    (Volts) 

2 

v'    e'    (Millivolts) 

0.0000 

2.780 

2.95 

12.4 

0.0625 

2.730 

2.95 

13.4 

0.1250 

2.630 

2.95 

15.6 

0.1875 

2.500 

2.94 

19.0 

0.2500 

2.330 

2.92 

22.6 

0.3125 

2.150 

2.90 

26.2 

0.3750 

1.960 

2.89 

29.4 

0.4375 

1.650 

2.85 

34.2 

Radial 
Position 
fin) 
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Table  8 
Thermocouple  Data  for  Re=30,700  Run 

0.00   0.05   0.10   0.15   0.20   0.25   0.30   0.35   0.40   0.45  NOISE 


/  e'2 
(Milli- 
volts)  20.0   21.0   23.5   28.5   30.5   34.5   38.5   42.5   45.0  47.5   8.0 


Table  9 

Thermocouple  Data  for  Re=41,500  Run 

Radial 
Position 

(in)    0.00   0.05   0.10   0.15   0.20   0.25   0.30   0.35   0.40   0.45  NOISE 


/e'2 
(Milli 
volts)   15.0   15.7   18.3   21.5   25.0   28.5   32.0   35.5   38.0   39.5   6.0 


Table  10 
Thermocouple  Data  for  Re=49,800  Run 

Radial 

Pos  it  ion 

(in)    0.00  0.05   0.10   0.15   0.20   0.25   0.30   0.35   0.40   0.45  NOISE 


v  e'^ 
(Milli- 
volts)  12.0   12.8   14.8   17.5   20.5   24.0   27.0   30.0   32.0   33.2   3.0 
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Table    11 
Non- Isothermal   Hot-WireMBef  ore"Calibrat  ion   Data 


AP  (in)- 


R    -R    (Ohm) 
w      a       ± 

2.20 

1.620 

1.330 

0.932 

0.480 

1.00 

1.865 

1.820 

1.795 

1.740 

1.665 

1.20 

1.990 

1.940 

1.915 

1.850 

1.775 

1.40 

2.125 

2.070 

2.045 

1.955 

1.875 

1.60 

2.230 

2.175 

2.145 

2.080 

1.960 

1.80 

2.325 

2.265 

2.235 

2.180 

2.075 

2.00 

2.410 

2.355 

2.320 

2.265 

2.150 

2.20 

2.500 

2.435 

2.400 

2.340 

2.225 

2.40 

2.575 

2.515 

2.475 

2.410 

2.295 

2.60 

2.650 

2.590 

2.550 

2.490 

2.360 

2.80 

2.725 

2.660 

2.620 

2.555 

2.425 

3.00 

2.790 

2.725 

2.680 

2.620 

2.485 

3.20 

2.855 

2.785 

2.745 

2.686 

2.585 

3.4  0 

2.920 

2.84  5 

2.805 

2.746 

2.640 

'-'•'  Orifice  Meter  Manometer  Fluid  SPG  2.95 
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Table    12 
Cold  Resistances    in" 3ef ore"Calibrat ion 

P    (in.)*  R      uCOHM) 

cold 


2.120 
1.620 
1.330 
0.932 
0.480 


7.42 
7.46 
7.49 
7.57 
7.80 
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Table    13 
Temperature   Difference   Between  Wall    and    Centerline 

AP    <tn>*  Tr/R=l    -    Tr/R  =  0      ^ 
2.120  46 

1.620  51 

1.330  54 

.932  61 

.480  74 

-Orifice  Meter  Manometer  Fluid  SPG  2.95 
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Table     14 

N on- Isothermal   Hot-Wire    RMS    Voltage 
Fluctuation   Data    for     r/R  =  0 


AP  (in)" 


R    -i 
w 

(Ohm) 
+ 

1 

.00 

1 

.20 

1 

.40 

1 

.60 

1 

.80 

2 

.00 

2 

.20 

2 

40 

2 

60 

2 

80 

3 

00 

3. 

20 

3. 

40 

2.120 


24.8      Milli- 
volts 
24.8 

24.7 

24.5 

24.3 

24.3 

24.2 

24.1 

24.1 

24.0 

24.0 

23.9 

23.9 


1.620      1.330      0.932      0.480 


27.0  28.8  31.6  38.3 

26.9  28.4  31.0  37.2 

26.3  28.2  30.3  36.2 
26.0  27.9  29.8  35.3 
25.7  27.4  29.2  34.3 
25.5  26.9  28.9  33.6 
25.2  26.7  28.5  32.9 
25.0  26.4  28.1  32.0 
24.9  26.2  27.7  31.4 
24.7  25.9  27.5  31.0 
24.5  25.8  27.2  30.6 
24.5  25.6  27.0  30.0 

24.4  25.5  26.7  29.8 


*  Orifice  Meter  Manometer  Fluid  SPG  2.95 


Table     15 

Non- Isothermal    Hot-Wire   RMS   Voltage 

Fluctuation  Data  for  r/R=.5 
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AP(ln)*     > 

R    -R    (Ohm) 
w      a 

2.120 

1.620 

1.330 

0.932 

0.480 

Milli- 

1.00 

42.9 

volts 

45.0 

46.8 

50.5 

59.7 

1.20 

42.9 

44.8 

46.8 

50.0 

58.3 

1.40 

42.6 

44.6 

46.2 

49.2 

57.1 

1.60 

42.6 

44.0 

45.8 

48.7 

56.2 

1.80 

42.7 

44.0 

45.6 

48.0 

55.0 

2.00 

42.7 

43.8 

45.1 

47.5 

54.2 

2.20 

42.3 

43.5 

44.9 

47.0 

52.8 

2.40 

42.3 

43.3 

44.8 

46.2 

52.2 

2.60 

42.2 

43.2 

44.5 

46.0 

51.7 

2.80 

42.8 

43.2 

44.5 

45.8 

50.8 

3.00 

42.6 

43.0 

44.3 

45.2 

50.3 

3.20 

42.6 

43.0 

44.3 

45.1 

49.5 

3.40 

42.6 

43.0 

44.0 

44.8 

4  9.2 

Orifice  Meter  Manometer  Fluid  SPG  2.95 
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Table     16 
Non-Isothermal   Hot-Wire"Af ter"Calibrat ion  Data 


AP 


(in)* 


R    -R    (Ohm) 
w      a 

2.120 

1.620 

1.330 

0.932 

0.480 

1.00 

2.010    (Volts) 

1.94  5 

1.920 

1.860 

1.775 

1.20 

2.140 

2.100 

2.070 

1.985 

1.890 

1.40 

2.235 

2.215 

2.180 

2.120 

2.000 

1.60 

2.365 

2.320 

2.285 

2.220 

2.125 

1.80 

2.470 

2.420 

2.385 

2.320 

2.215 

2.00 

2.565 

2.515 

2.480 

2.405 

2.295 

2.20 

2.655 

2.600 

2.565 

2.495 

2.375 

2.40 

2.740 

2.685 

2.645 

2.570 

2.450 

2.60 

2.820 

2.765 

2.720 

2.645 

2.525 

2.80 

2.900 

2.840 

2.795 

2.720 

2.595 

3.00 

2.970 

2.915 

2.860 

2.785 

2.655 

3.20 

3.070 

2.980 

2.930 

2.845 

2.715 

3.40 

3.135 

3.075 

2.995 

2.910 

2.775 

'-'•'  Orifice  Meter  Manometer  Fluid  SPG  2.95 


56 


Table 

17 

Cold   Resistances 

in* 

'Aft 

er"Calibration 

P    (in.)" 

.•, 

R         .     (OHM) 
cold 

2.120 

7.21 

1.620 

7.24 

1.330 

7.26 

0.932 

7.29 

0.480 

7.43 

APPENDIX    B 
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APPENDIX  B 
ORIFICE  l^ETER  CALIBRATION 

The  orifice  meter  used  was  calibrated  in  place  against  a  pitot  tube 
located  at  the  centerline  of  the  test  section.   von  Karman's  universal 

1  Q 

velocity  distribution   was  used  to  calculate  flow  rates  from  the  pitot 
tube  velocity  measurements.   The  calibration  curve  is  shown  in  Figure  14 

It  is  found  as  anticipated  that  a  plot  of  flow  rate  versus  the  square 

19 
root  or  the  pressure  drop  is  a  straight  lint. 
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Figure  14.   Orifice  Meter  Calibration  Curve 
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APPENDIX  C 

AMPLIFIER  FREQUENCY  RESPONSE 
The  frequency  response  of  the  Tektronix  amplifier  is  shown  on 
the  following  page. 
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APPENDIX   D 
RESISTANCE-TEMPERATURE    CURVE    FOR    THE    HOT-WIRE 

Figure    16    shows    the    relationship   between   resistance    and    temperature 

for    the    hot    wire    used    in    the    non- isothermal    flow   measurements.      The 

3 
dashed    curve    is    the    approximate    relationship    suggested    by    Kronauer . 
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APPENDIX   E 

CALIBRATION    OF    HOT-WIRE    IN    ISOTHERMAL  FLOW 

The    isothermal    flow   data   were    taken    on    two   different    days.      The    set 
of    data    for    a   Reynolds    Number    of   49,800  was    taken    on    the    first    day.      The 

calibration    curve    for    these    data    is    shown    in   Figure    17.      The    constant    C 

-2 
in   Equation    18,    which    is    the    slope    of    the    line    in   Figure    17    is    1.67x10 

The    data    at    Reynolds    Numbers    of    30,700   and    41,500  were    taken    on    the 

second    day.      The    calibration    curve    for    these    data    is    shown    in   Figure    18. 

-2 
The    constant    C    for    these    data    is    1.69x10 

-2 
The    value    of   C    calculated    from  Equation    17    is    0.84    x    10       .      The 

discrepancy    between    this    value    and    that     obtained    from   calibration    is    due 

to  wire   history   effects    in     a    and   R      and    the    error    inherent    in    King's 

o 

potential    flow    value    for    B. 
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Figure    17.       Isothermal   Calibration    of   Hot-Wire    I 
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Figure    18.       Isothermal   Calibration    of   Hot-Wire    II 
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APPENDIX  F 

LISTING  OF  CURVE  FITTING  PROGRAM 
The  computer  subroutine  used  to  perform  least  square  curve  fitting 
calculations  is  listed  on  the  following  pages. 
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APPENDIX  G 

CALIBRATION  OF  HOT-WIRE  IN  NON- ISOTHERMAL  FLOW 

Figures  19  through  23  represent  the  calibration  data  used  to  compute 

the  results  presented  for  non- isothermal  flow.   The  curves  drawn  through 

the  data  represent  the  least  squares  fit  of  parabolas  to  the  data.   These 

parabolas  were  differentiated  analytically  to  obtain  the  values  of  S   and 

u 

S   presented  in  Figures  9  and  10. 

Figures  24  through  26  show  calibration  curves  for  representative 

values  of  u  and  T   -  I   before  and  after  fluctuation  measurements  were 
w     a 

made.   Table  18  gives  a  complete  tabulation  of  S   and  S   as  functions  of 

u  and  T   -  T   before  and  after  the  fluctuation  measurements, 
w    a 
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Figure    19.      Non-Isothermal   Calibration    of   Hot -Wire    I 
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Figure  20.   N on- Isothermal  Calibration  of  Hot -Wire  II 
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Figure   21.      Non- Isothermal   Calibration   of   Hot-Wire    III 
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Figure    22.      Non-Isothermal  Calibration    of   Hot-Wire    IV 
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Figure  23.   Non-Isothermal  CalibraLion  of  Hot-Wire  V 
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Figure    24.      Calibration   Drift    I 
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Figure  25.   Calibration  Drift  II 
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Figure  26.   Calibration  Drift  III 
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APPENDIX  H 

KOVASZNAY  FLUCTUATION  DIAGRAMS 

Figures  27  through  29  show  the  Kovasznay  Fluctuation  Diagrams  for 
three  representative  cases. 

The  parabolas  that  fit  these  data  are: 


F  igure 
27 

28 

29 


Parabola 


,2   — -  S  S   2 

(  ^~j   -    t,Z)  =  2.64  +  16. 9(^)  +  19. M^) 
S„,  T         bT 


,2 


S   2 


—  ■  t'2)  =  0.43  +18. 9^)  +  16. 2(-^) 

s/  bT  bT 


,2  S  S   2 

(  ^-J  -    t,Z)  =  0.21  +  42. 6^)  +  47. 3^) 
S  ^  bT  bT 


u  t    is  then  found  by  dividing  the  coefficient  of 


u 


by  -2, 


,2 


T  S   2 
is  the  square  root  of  the  coefficient  of   (   "z- )  .   The 

bT 


and  v    u 

results  from  the  parabolas  given  above  are  listed  in  Table  19 


Table  19 


,2 


Sample  Results  for  *  u   and  u  t 


Reynolds 

Mumber 

J      - 2 
v    u 

And    r/R 

u't' 
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0.0 

4.36 

-8.45 

38,200 

0.0 
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-9.45 

38,200 
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Figure    27.      Kovasznay   Fluctuation   Diagram   I 
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Figure   28.      Kovasznay  Fluctuation  Diagram   II 
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Figure  29.   Kovasznay  Fluctuation  Diagram  III 
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APPENDIX  I 
ERROR  ANALYSIS 

The  results  of  the  isothermal  flow  measurements  with  the  hot-wire 
deviate  from  the  results  of  Sandborn  "  and  Laufer   by  less  than  10 
percent . 

The  results  of  a  preliminary  non- isothermal  run  duplicating  the 
conditions  of  Rodriguez-Ramirez  and  using  the  rapid  response  thermocouple 
agreed  with  his  results  within  5  percent . 

Since  the  errors  in  the  results  of  the  non- isothermal  flow  measure- 
ments with  the  hot-wire  are  quite  large,  a  detailed  estimate  of  their 
magnitude  is  needed. 


that  the  errors  in   V  u '  /u 


It  is  assumed  that  the  errors  in   \u'  /ut  and  u't'/u   (T  -T, ) 

t  tot 

are  primarily  the  result  of  errors  in  S   and  S_ .   In  order  to  obtain  an 

u       T 

estimate  of  the  errors  in  S   and  Sm ,  an  analysis  was  made  of  the  dis- 

u      I 

crepancies  between  successive  calibrations  as  tabulated  in  Table  18. 

Note  that  the  actual  errors  in  S   and  S_  come  from  a  continuous 

u      T 

drift  in  E  during  a  given  calibration  as  outlined  in  the  DISCUSSION 
section.   Since  it  is  not  possible  to  measure  this  drift,  recourse  was 
made  to  analyzing  by  this  method.   Thus  the  error  estimates  obtained 
below  do  not  represent .confidence  limits  on  the  results  but  only  rough 
approximations  to  the  confidence  limits.   In  Table  18  ,  the  largest 


86 


discrepancies  occur  for  u  =  137  ft/sec.   The  percentage  error,  based  on 

the  before  value,  between  the  before  and  after  values,  was  determined 

for  each  entry  in  the  137  ft/sec  column.   All  the  S   errors  were  summed 

and  divided  by  the  number  of  entries  and  likewise  for  S  .   Thus,  average 

u 

errors  were  determined  to  be  19  percent  in  S   and  14  percent  in  S  . 

The  propagation  of  these  errors  into  the  coordinates  of  the  Kovasznay 
Fluctuation  Diagram  is  as  follows: 


2     2 
A  (  £2  -   t  ) 

ST =    2  AST     =  28  percent 


?      2  ST 

e'  •  -  t' 


A(Su) 


ST    =ASu   +  AST   =  33  percent 


s      s      s 

u         u        T 


ST 


These  errors  were  assumed  to  apply  to  all  values  of  u  and  r/R. 

S 
The  value  of  —   for  each  point  of  the  Kovasznay  Fluctuation 

bT 
Diagram  for  u  =  66  ft/sec  r/R  =  0  was  increased  by  33  percent,  and  the 


2      2 
value  of  e '   -  t'    was  decreased  by  28  percent.   A  parabola  was  fitted 

*T2 


to  the  resulting  points  and  u'    and  u't1   determined.   The  new  values 
differed  from  the  original  values  by  72  percent  and  100  percent  respec- 
tively.  Since 


Vu.2  =  i .  u. 


i   A 


"7?      2  -   — 

u  u 


the  error  in  \  u '    i 


s  36  percent.   This  36  percent  error  was  used  as 


the  limit  of  error  for  all  values  of  \u'  /u   ,  and  100  percent  was  used 


V  -1 


for  all  values  of   u't'/u  T (T  -T  ) 
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